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We present a quantum version of the vacuum Bianchi IX model by implementing affine coherent 
state quantization combined with a Born-Oppenheimer-like adiabatic approximation. The analytical 
treatment is carried out on both quantum and semiclassical levels. The resolution of the classical 
singularity occurs by means of a repulsive potential generated by our quantization procedure. The 
quantization of the oscillatory degrees of freedom produces a radiation energy density term in the 
semiclassical constraint equation. The Friedmann-like lowest energy eigenstates of the system are 
found to be dynamically stable. 


Introduction. The Friedmann-Robertson-Walker model 
is successfully used to describe the data of observational 
cosmology (see e.g. BE])- Nevertheless, the isotropy 
of space is dynamically unstable towards the big-bang 
singularity [3]. On the other hand, if the present Uni¬ 
verse originated from an inflationary phase, then the pre¬ 
inflationary universe is supposed to have been both in¬ 
homogeneous and anisotropic. As evidence suggests (see 
m, 0), the dynamics of such universe backwards in time 
becomes ultralocal and effectively identical with the ho¬ 
mogeneous but anisotropic one at each spatial point. In 
both cases quantization of the isotropic models alone ap¬ 
pears to be insufficient. Hence the quantum version of 
an anisotropic model, comprising the Friedmann model 
as a particular case, is expected to be better suited for 
describing the earliest Universe. 

In the Letter we advocate a new quantization method 
of the dynamics of a vacuum Bianchi type IX geometry, 
the Mixmaster universe. We identify a soluble sector of 
this model, which lies deeply in the quantum domain and, 
as we show, contains relevant physics. 

The Mixmaster universe exhibits a complex behavior 
[5]. As it collapses, the universe enters chaotic oscilla¬ 
tions producing an infinite sequence of distortions from 
its spherical shape [7]. Those distortions essentially cor¬ 
respond to the level of anisotropy and may be viewed as 
an effect of a gravitational wave evolving in an isotropic 
background [B] . Dynamics of this wave is nonlinear, and 
its interaction with the isotropic background fuels the 
gravitational contraction. Not surprisingly, the quan¬ 
tization of Bianchi IX model is a difficult task. Some 
formulations can be found in literature, including the 
Wheeler-DeWitt equation [5] or, more recently, a formu¬ 
lation based on loop quantum cosmology imn]. How¬ 
ever, the search for solutions within these formulations 
is quite challenging m [12] leaving the near big-bang 
dynamics largely unexplored. The most recent develop¬ 
ments, e.g. m, do not address the singularity resolution. 

To fill this gap we propose a quantum Mixmaster dy¬ 
namics, which originates from affine coherent state (ACS) 


quantization that was recently used to obtain the quan¬ 
tum Friedmann model m. It was shown that ACS 
quantization causes some new terms to appear in the 
quantum Hamiltonian, producing a strong repulsive force 
counteracting the contraction of universe. The capacity 
to resolve the singularity constitutes the basic advantage 
of our quantization method. In order to solve the dy¬ 
namics in the present, more complex setting, we employ 
the adiabatic approximation widely utilized in quantum 
molecular physics [miis]. This approach is reasonable 
when the vibrations of the shape of the universe are sig¬ 
nificantly faster than the contraction of its volume. 

The main result is a semiclassical Friedmann-like equa¬ 
tion obtained from expectation values in ACS, a descrip¬ 
tion peculiar to our approach. In that equation, the 
expansion of the universe is governed by two terms of 
quantum origin. The first one is proper to the quantum 
Mixmaster model and corresponds to the energy of the 
wave in an eigenstate. It is proportional to the energy 
level number. The other one, which is more universal, 
corresponds to the repulsive potential preventing the sin¬ 
gularity. The lowest energy eigenstates of this system are 
interpreted as the quantum Friedmann universe supple¬ 
mented with vacuum fluctuations of the anisotropy. 

Beyond issues of singularity resolution, the Friedmann- 
like equation describes two novel and rather surpris¬ 
ing properties of the quantum dynamics. Firstly, the 
anisotropic degrees of freedom remain in their lowest en¬ 
ergy states during the quantum phase consistent with our 
approximation. This implies that the quantum Fried¬ 
mann model, unlike its classical counterpart, is in fact 
stable with respect to the anisotropy. Therefore, the clas¬ 
sical chaos is suppressed within the considered domain. 
Secondly, during the contraction the quantum energy of 
anisotropy grows much slower than it does on the clas¬ 
sical level. Namely, it effectively gravitates as radiation 
leading to a significant reduction in the overall gravita¬ 
tional pull from anisotropy due to quantum effects. 

Classical Hamiltonian. The Friedmann equation ex- 
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tended to anisotropic universes, with c = 1 = SttG, reads: 




r Q „ rr 

- R- = -p, 


( 1 ) 


where H = a/a (dot denotes the derivative with respect 
to the cosmological time) is the expansion parameter, p 
is the energy density of matter, and is the spatial 
curvature. The additional term is the total shear of 
the spatial section and is non-vanishing for anisotropic 
models. Due to its negative sign, the shear drives the 
gravitational collapse and it eventually dominates the 
dynamics. For this reason we neglect the attraction of 
matter by putting p = 0, that is, we restrict our consid¬ 
erations to the vacuum Bianchi IX, whose dynamics near 
the singularity is effectively the same as in the presence 
of perfect fluid [1]. 

The Mixmaster describes the space-time metric ds^ = 
—dt^ -I- jcr'-a^, where a is an averaged scale factor 

and cr* are differential forms on a three-sphere (covering 
group of the rotation group) satisfying dcr* = ^Cijkcr^ AaR 
The diagonal form of the metric is assumed '■= 

diag 6"*^+), where j3± are dis¬ 

tortions parameters [6]. 

In terms of these variables, the shear is the kinetic 
energy of anisotropic distortion: T? = (p^ -|-p^)/24a®, 
where the momenta p± are canonical conjugates to /3±. 
The spatial curvature ^R grows due to the overall con¬ 
traction of space, but decreases due to the growth of 
anisotropy. This last circumstance leads to the back- 
reaction from spatial curvature on the shear and, as a 
result, oscillations in I3± occur. As there is no mat¬ 
ter content in our model, /3± describe a sort of gravi¬ 
tational wave. The curvature can be split into isotropic 
and anisotropic parts: ^R = 3(1 — V{fi))/2a?, where 
F(/3) is the anisotropy curvature potential 0: 


F(/3) = (^{e-^^+ - 2cosh(2y3/3_))^ - 4^ -b 1. 

As shown in Fig. this potential has three “open” C 3 „ 
symmetry directions. One can view them as three deep 
“canyons”, increasingly narrow until their respective wall 
edges close up at the infinity whereas their respective bot¬ 
toms tend to zero. Due to its (almost) confining shape, 
V is expected to produce a discrete energy spectrum on 
the quantum level. 

The generalized Friedmann equation Q may be 
rewritten as 




1 1^2 V(/3) 

4a^ 6 4a^ 


( 2 ) 


where the isotropic background geometry on the l.h.s. 
is pulled by the energy of anisotropic oscillations. The 
energy of oscillations scales roughly as a“®. 

It follows that the Hamiltonian constraint to be quan¬ 
tized reads in canonical variables: 




( 3 ) 



FIG. 1. Global picture of the potential V(/3) near its mini¬ 
mum. Boundedness from below, confining aspects, and three 
canyons are illustrated. 


where q = and = I6d^a are more suitable to ACS 
quantization, and where 

R, = ^{pl+pl) + \q^'^V{p) (4) 

is the anisotropy energy. The closed Friedmann- 
Robertson-Walker (FRW) geometry is obtained by 
putting p± = 0 and j3± = 0, or simply Hq = 0. 

The Hamiltonian constraint ([^ resembles a diatomic 
molecular Hamiltonian with the pairs {q,p) and (/3±,p±) 
playing the role of the reduced nuclear and electronic 
variables, respectively. In molecules, the motion of nuclei 
is slow enough in comparison with electrons so the motion 
of electrons may be approximated as becoming immedi¬ 
ately adjusted to varying positions of nuclei. However, 
the coupling between the nucleus-like and electron-like 
degrees of freedom in the present model differs from the 
usual molecular case for which the validity of the approx¬ 
imation rests upon the smallness of the ratio between the 
nuclei and electron masses. In the present case, described 
by Eqs. ([^ and Q, the “mass” of the degrees of freedom 
/3± evolves as q. Thus, it goes to zero near the singu¬ 
larity, (7 = 0. On the other hand, the “mass” of the 
degree of freedom q in Eq. ^ is constant. Thus, close 
to singularity the latter may be regarded as “heavy” in 
comparison with the anisotropic variables that can be 
treated as “light”. 

Quantum Hamiltonian. The six-dimensional phase space 
of the Mixmaster universe is quantized as follows: (A) 
The isotropic variables form the canonical pair {q,p) liv¬ 
ing in a half-plane. That half-plane can be viewed as 
the affine group. We resort to one of its two unitary ir¬ 
reducible representations, denoted by U, to build from 
a suitable fiducial vector \v) (where > 0 is a free pa¬ 
rameter) a family of affine coherent states (i.e., wavelets) 
\q,p) := U[q,p)\y). These ACS’s are then used to consis¬ 
tently quantize the isotropic variables. While the method 
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provides the usual p = —ihdq, and q defined as the multi¬ 
plication by g, its interest lies in the regularization of the 
Hamiltonian m- This approach together with \v) was 
introduced for cosmological models in M- Next, we use 
the ACS’s to obtain a semiclassical description, which 
enables us to analyze the effective dynamics of isotropic 
variables. (B) For the anisotropic variables, each canon¬ 
ical pair (/3±,p±) lives in the plane. Thus, it is natural 
to proceed with the usual canonical quantization which 
yields p± = —ihdp^ , and j3± being the multiplication by 

P±- 

The quantized Hamiltonian corresponding to ([^ and 
issued from quantizations (A) and (B) above reads 

where 

+ ( 6 ) 

The A^’s are purely positive numerical constants depen¬ 
dent on the choice of the ACS. With the choice made in 
our previous paper m all these constants are simple ra¬ 
tional functions of modified Bessel functions Ki{v). We 
note in ([^ the appearance of the repulsive centrifugal 
potential term It is the signature of the ACS 

quantization, which is consistent with the half-plane ge¬ 
ometry, and it regularizes the quantum Hamiltonian for 
small q m- As the universe approaches the singularity, 
g —>■ 0, this centrifugal term sharply grows in dynamical 
significance. 

We consider the oscillations of /3± fast in comparison 
with the contraction of the universe. It legitimates the 
adiabatic approximation, in a way analogous to the Born- 
Oppenheimer approximation (BOA) |15[II6) widely used 
in molecular physics. Due to the confining property of 
V, the operator Hg at fixed q has a discrete spectrum. 
In accordance with BOA we assume that the anisotropy 
degrees of freedom (3± are frozen in some eigenstate of 
Hq with eigenenergy (A^ = 0,1,...) evolving adia- 
batically. Thus, the light degrees of freedom (3± can be 
averaged leading to the Hamiltonian: 

Focussing on the deep quantum domain, we look at the 
first energy levels near the ground state of Hq. Therefore 
we essentially investigate the domain near the minimum 
of V(P). Within the harmonic approximation V{/3) « 
8(/3+ -f 13^) near its minimum /3± = 0, the eigenener- 

gies are found to be Cq^'^ ~ tiq~'^/^\/ 2 ^ 2 ^z{N + 1) with 
N = n_|_ + n_. The quantum numbers n± correspond 
to the independent harmonic oscillations in j3+ and /3_. 
Consequently, the approximation of the eigenvalues of T-Lq 


may be written as |18j : 

eW~^v/2^(iV + I). (8) 

The expression for is rather a rough approximation 
for large values of N, since C(/3) is highly nonharmonic 
far away from its minimum. But for small values of N^ 
this expression is valid at any value of q. 

Validity of approximation. We notice in Eq. ([^ that the 
discrete spectrum part in Eq. Q becomes a small pertur¬ 
bation at large g, a range for which BOA possibly loses 
its validity, whereas it gains all its value at small g. From 
the mass criterion, our approach based on BOA is legit¬ 
imate as g assumes its values near the singularity g = 0. 
Furthermore, our procedure of quantization generates a 
supplementary repulsive potential that prohibits the sys¬ 
tem to access the singularity neighborhood g S {0,qm) 
with some very small bound qm > 0, which depends on 
the initial state. 

Furthermore, calculations made in molecular physics 
beyond BOA (the so-called vibronic approximation) show 
that the mass criterion is in fact too strong: a significant 
breakdown of BOA only occurs when different eigenen¬ 
ergy curves g i—>■ Cq^'^ of Hq are crossing. In our approach 
these crossings do not occur, at least for the lowest levels 
of Hq in Eq. (|^. 

This reasoning based on molecular physics is robust, 
but qualitative in our case, due to the coupling between 
the g and degrees of freedom which is not of molec¬ 
ular type. In [18] we weaken the adiabatic condition by 
allowing the quantized oscillations to be excited by the 
semiclassical dynamics of the isotropic background de¬ 
scribed below, for a fixed N. We find that the excitation 
is indeed very limited, which justifies our approach. 
Semiclassical Hamiltonian. We introduce a semiclassical 
observable associated with the quantum Hamiltonian 0 
as its expectation value CA{q,p) ■= {q,p\CA\Q,p) in the 
ACS state \q,p) peaked on the classical phase space point 
(g,p) in the half-plane. 


CAiq,p) = ^ (p' + + 1) , 

where the ’s are positive numerical constants [I8j which 
are also simple rational functions of modified Bessel func¬ 
tions With our choice of |i^), at large v, ^ 1, 

* ^ 4 and ^ vj^. For the consistency of our semi¬ 

classical description we have rescaled the fiducial vector 
so that (g,p|g|g,p) = g and (g,p|p|g,p) =p. 

The Hamiltonian constraint imposed at the semiclas¬ 
sical level, CA[q,p) = 0, leads to the semiclassical 
Friedmann-like equation: 






a 


■^5 ^ c W SttGH 

® 4 Va/ 3c 


(7V + I)J. (9) 
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FIG. 2. Three periodic semiclassical trajectories in the half¬ 
plane (a, H) from Eq. They are smooth plane curves. We 
use standard units SttG = c = = 1 and choose i/ = 1. Blue 

dotted curve for = 0, green dotdashed for = 1 and red 
dashed for A'" = 2. 


where we have restored physical constants and the stan¬ 
dard cosmological variables. 

The above semi-classical constraint admits smooth tra¬ 
jectories for all values of N only if ly G (0,7.19). For 
ly > 7.19, Eq. (|^ has no solution for the smallest values 
of N. The solution of (|^ for a is a periodic function, 
a G [o_,a+] with a_ > 0 and a+ < oo, and resolves the 
cosmological singularity of the Mixmaster universe. In 
Fig. 2 we plot a few trajectories in the half-plane (o, H). 
The classical closed FRW model is recovered at h = 0 
and large values of ly. 

Discussion. Our semi-classical analysis of the Mixmaster 
universe leads to the modified Friedmann equation (§. 
The left-hand side describes the isotropic part of geom¬ 
etry. The Hubble parameter squared is accompanied by 
the repulsive potential of purely quantum origin, which 
grows as during the contraction. At small volumes, it 
efficiently counteracts the attraction of common forms of 
matter, forcing the collapsing universe to rebound. The 
third term is the usual isotropic spatial curvature. 

The right-hand side of Eq. (§ describes the quantized 
energy of the anisotropy oscillations. The energy is dis¬ 
crete and increases linearly with integer Af, as expected 
in our harmonic approximation. Within the adiabatic 
approximation, that quantum number is conserved dur¬ 
ing the evolution. The energy of the quantum oscilla¬ 
tor evolves due to the a-dependent coefficients in front 
of its kinetic and potential terms given in Eq. ([^. The 
ratio between the coefficients determines the oscillator’s 
frequency, which is proportional to a^. The energy of 
the oscillations at the quantum level is multiplied by the 
frequency and consequently it scales as a~‘^. (This be¬ 
comes a poor approximation for high values of N, due 
to the breakdown of the harmonic approximation). It is 


quite a contrast to the classical wave, whose total energy 
is approximately unaffected by its time-dependent fre¬ 
quency, and therefore scales as a“®. Thus, the growth of 
the attractive force induced by the anisotropy is signifi¬ 
cantly reduced in the semiclassical dynamics. This essen¬ 
tial dissimilarity between the classical and semiclassical 
dynamics is due to the fact that on the quantum level 
the contraction of space is driven by a quantum average. 

Let us note that the energy of the wave does not vanish 
even in the ground state A^ = 0 due to the zero-point 
quantum fluctuations corresponding to the classical state 
/3± = p± = 0. In [T5] we go beyond the adiabatic ap¬ 
proximation to check if there is a significant excitation of 
the wave’s energy level during the semiclassical evolution 
of the background geometry. The method is essentially 
the same as the one used to discuss the generation of 
primordial power spectra in inflationary cosmology. We 
show that the wave in fact remains in its lowest energy 
states during the quantum phase. It confirms that the 
quantum FRW universe, unlike its classical version, is 
dynamically stable with respect to the small isotropy 
perturbation. Therefore, it seems that the quantum 
closed Friedmann model may be successfully used to 
describe also the earliest Universe, provided that the 
corresponding Hamiltonian is supplemented with the ef¬ 
fect of the zero-point energy generated by the quantized 
anisotropy degrees of freedom of the Mixmaster universe. 
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